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Testing fixed points in the 2D 0(3) non-linear a—model. 
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Using high statistic numerical results we investigate the properties of the 0(3) non-linear 2D 
n-model. Our main concern is the detection of an hypothetical Kosterlitz-Thouless-like ( I\T ) phase 
transition which would contradict the asymptotic freedom scenario. Our results do not support such 
a KT- like phase transition. 
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I. INTRODUCTION 


The bidimensional 0(N ) non-linear sigma model is defined on the lattice by the action 

S = -P X ■ <j>(x + n) . 
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Here cj)(x) is a iV-component vector constrained by cj) 2 (x) = 1. Perturbation theory predicts that this model is 
asymptotically free if TV > 3 |Q|^] . Large distance observables must scale accordingly with the renormalization group 
equation. In particular we find for the correlation length 
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and for the magnetic susceptibility 
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The coefficients ak and bk can be evaluated perturbatively starting from 3-loop order. They describe non-universal 
(scheme dependent) corrections to asymptotic scaling. The constants Cj and C x cannot be evaluated perturbatively. 
There are many possible definitions of correlation length. All of them scale in the same way, Eq. (§ with the same 
coefficients a*,. The constant Cf is instead definition dependent. The exponential fall-off of the correlation function 
at large distances provides the so-called exponential definition whose constant C% can be obtained exactly (alas not 
rigorously!) with the Bethe-Ansatz technique 
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C x is known at second order in the 1/N expansion g. In the ratio x/£ 2 the exponential dependence on (3 disappears, 
and the quantity 
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should stay constant for large enough (3. The coefficients Cfc can be easily calculated from ak and bk- 

According to the Mermin-Wagner theorem Q the continuous O(N) simmetry cannot be broken in two dimensions. 
Perturbative expansions are constructed around a non-symmetric trivial vacuum, so their validity is not guaranteed. 
By analysing the percolation properties of the Fortuin-Kasteleyn clusters [Q it has been argued |^,|| that perturbative 
methods are inadequate to describe the model. The main conclusion from these analyses is that the non-linear sigma 
model in 2 dimensions has no mass-gap, contradicting the results of and that for any Af > 2 it undergoes a 

KT like phase transition at a finite inverse temperature (3kt- 
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If this is the case, then the ratio 


Rkt = - P) r ( 6 ) 

should be asymptotically constant for (3kt — /3 positive and small. Here is a critical exponent. For the 0(2) model, 
where the KT transition is generally trusted, 77 = 1/4. On the other hand, renormalization group considerations and 
numerical evaluations of r yield the bound |r| < 0.1 (h] [l^] . 

It is well known that in the 0(3) model there is not a good numerical evidence of the onset of the asymptotic scaling 
regime. In juj the authors find that in this model the numerical data for the ratio Rkt with rj = 1/4 can be well 
fitted with a constant, while the ratio Rpt is clearly decreasing. They make use of data obtained from a simulation 
with the standard action Eq. ©■ 

In ||l 6 || it has been shown however that Rkt is not constant. To do this, the O-loop Symanzik improved action was 
numerically simulated at high statistics (O(10 7 ) measurements), including corrections to asymptotic scaling up to 3 
loops. In |uj we claimed that the discrepancy on the constancy of Rkt between |l(| and |L5| originates from the 
higher resolution of the data in ]n| . Arguments against these conclusions has been raised based on possible distorting 
effects of the antiferromagnetic coupling present in the improved action. Here we want to perform an analogous study, 
at high statistics too, on the 0(3) model with standard action where such a coupling is absent. 

For completeness we will show the results for Oj, C x and the Rpt ratio obtained from our data. 


II. MONTE CARLO SIMULATION 

We have simulated the standard action Eq. (ilj) on lattice sizes L by using the Wolff algorithm (jT^. The wall-wall 
correlation function 


G(i) = 7 x )) 
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has been measured with improved estimators (JT 8 J . We have always chosen L big enough to avoid sizable finite volume 
effects: L/£ > 7 

Using these data we can evaluate the correlation function at momentum p 


G(p) = Lj2e ipt G(t), 


( 8 ) 


the magnetic susceptibility x = G(0) and the correlation length from the exponential fall-off of Eq. 

£= “/i™ In G(t)' 

The data for these quantities are listed in Table [j], together with the energy density defined as 

E = (0(0) • 0(0 + A)) 

where A is a generic direction and there is no a sum over that p. 
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III. RESULTS 

Although the present note is mainly concerned on the existence of a phase transition at a finite temperature, we 
give for completeness also the results for 6 \, C x and Rpt- 

By using Eq. (!2j) and our Monte Carlo data we obtain a prediction about Cj, which we call CY c . In Fig. [I] we 
plot with empty symbols the ratio C^ ic /C$ versus /?, adding scaling corrections up to four loops | 2 (J (we have taken 
into account the corrections of |Hf). It is apparent that we are not in an asymptotic scaling regime, although the 
non-universal scaling corrections provide a systematic improvement. In the best (4-loop) case the lack of asymptotic 
scaling is around 15%. 

We have also used an effective scheme Q defined by the perturbative expansion of the energy Eq. (0) 
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The results in this scheme are displayed as full symbols in Fig. |T]. Here the lack of asymptotic scaling is still large: 
about 7-8% in the best case. 

Starting from Eq. (|j) we can study in the same way the ratio C X G /C x , using for C x the 0(1/N 2 ) result J|] 

C x /* V ^ = 0.0127. The percentage of discrepancy in Fig. [| between the Monte Carlo result and the (1/iV 2 ) analytical 
prediction is similar to that of the correlation length in Fig. [j]. It is interesting to note that the 3-loop correction is 
irrelevant, while the 4-loop one is not. 

The data in the effective energy-scheme for the Fig. || are displayed as full symbols. Here we see an excellent 
agreement between the Monte Carlo data at 4-loops and the 1/IV 2 approximation. We do not know how good is the 
1/JV 2 approximation to C x . Therefore we avoid drawing optimistic conclusions about the approach to the asymptotic 
scaling regime in this case. 

In Fig. we plot the logarithm of the ratio Rpr versus the logarithm of the correlation length. Again the full 
(open) symbols correspond to the effective (standard) scheme. We see that data cannot be considered constant, 
although corrective terms give a systematic improvement. We stress the fact that further corrections and further 
statistics do not worsen the results. From the larger (3 and considering the errors coming altogether from the 2-3-4 
loop corrections, we get In (C x /C^ 2 ) = 4.5705(39) which, by using Eq. ([|), leads to the prediction C x = 0.01507(6). 
This value however must be regarded as an upper bound because data in Fig. || are still descending. 

In Fig. [I] we plot the main result of the present paper: the ratio Rkt versus £ in a log-log scale. As it already 
happened for the Symanzik action, the shape of the curve indicates a non-constant ratio, contradicting the claims 
of [JHiJ. In calculating this figure we have omitted the power (/3 KT — P) r , see Eq. (||). Had the model to undergo a 
phase transition at finite coupling Pkt, the correlation length £ should behave like 


1 =■ 4exp (■7db?) - (12) 

If we try to fit our data to such functional form we obtain very unstable results (including more or less data points 
in the fit, the results for A, B and Pkt vary strongly). The fit to the form Eq. (|^) is much more stable. This 
phenomenon has a simple interpretation: Pkt is so large that Eq. © can be approximated to 

£ ss A' exp (B' (3) A 1 = A exp f E— \ B' = —(13) 

\VPktJ 2 (3]i-l 

and therefore we are fitting actually the combination B/P^t- If Pkt is so large and recalling that |r| < 0.1, then the 
power (Pkt — (3) r is almost constant within the narrow interval of our working p’s. This is why we omitted such a 
factor in the calculation of Rkt- 

Notice that the shape of the curve for Rkt can be obtained by using the perturbative expressions in Eqs. (fjj3). 
Owing to the incomplete asymptotic scaling, we obtain a 20% difference in the overall scale in the best case (4 loops 
and energy scheme). However the position of the minimum is well predicted: the perturbative expansion at 4 loops 
and energy scheme display a minimum at ln£ = 3.2 while the curve in Fig. 4 has minimum at ln£ = 3.5(4) (the error 
has been estimated by interpolating the top and bottom values of the error bars in the data of Fig. 4). 

The prediction of j||| is that the O(N) model behaves like the 0(2) model for N > 3. Therefore the set of critical 
exponents of 0(2 ) should also apply for 0(3). However one could imagine a set of slightly different exponents (after 
all, the 0(3) and 0(2) models have different numbers of degrees of freedom). Following this idea, we have made 
further fits on the data for £ to functional forms similar to Eq. ( 0 ) but with other powers of ((3kt — /?). In all cases 
the only sensible conclusion was that Pkt is very large. One can also try other versions of the ratio Rkt by varying 
7]. The tendency is that for larger (smaller) r] the curve acquires a positive (negative) slope for all points. We could 
not find a value of for which the whole set of points in Fig. [| flattens out. 


IV. CONCLUSIONS 

We have done a Monte Carlo simulation for the standard action of the 0(3) non-linear er-model in 2-dimensions. 
We have improved the statistics with respect to previous work taking advantage of the recently calculated corrections 
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to scaling |^o|,^lJ and energy |jn|. We have tested the perturbation theory predictions in both the standard and 
effective schemes. 

In Figs. 1-3 we have plotted our results for C)r, C x and Rpp as a function of /?. We see that the use of an effective 
scheme improves the asymptotic scaling in a sensible way. However, in the best case, we are still far from asymptotic 
scaling by roughly 10%. Much closer approaches to the scaling region were obtained by using the Symanzik action | {i(S| . 

It is remarkable that in all cases, the 4-loop correction in the effective scheme is negligible and that if we trust the 
0(l/iV 2 ) estimate of C x then our data have reached the scaling regime after including 3 or 4 loops. 

As for our main result, Fig. [I|, it does not give any support to the AT-scenario. The plot has the same shape than 
for the Symanzik action Jl6| | and this indicates that distorting effects due to the antiferromagnetic coupling are quite 
unlikely. In fact, in [fl(i| we worked with £ > 16 in order to avoid any strong coupling effect. 

As already shown in [^b|, the ratio Rpr is not constant but further corrections straighten it. This tendency is not 
satisfied for the Rkt ratio if we increase the statistics. In conclusion, we think that the main claim in |l^| cannot be 
maintained within the present status of numerical simulations because either (3kt is too large to see any effect of the 
A'T-transition or even it is infinity. 
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VI. NOTE ADDED 

To obtain the results displayed in Figs. 1-3, we have used the corrected values for the finite integrals published 
in ||^] needed for the 4-loop beta function. We notice that these corrected values do not change appreciably the 
results shown in 0- For instance, for the 0(8) model with standard action and energy scheme it was shown in 
that C^ c differs from Eq. (Q) by 0.5%; after the corrections of Q this difference becomes 0.4-0.5%. 
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FIG. 1. The ratio between non-perturbative constants C^ ic /C^ for the model. Empty circles (squares, diamonds) stand 
for the standard scheme 2-loop (3-loop,4 loop) results; filled circles (squares, diamonds) stand for the energy scheme 2-loop 
(3-loop,4-loop) results. 
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FIG. 2. 
of Fig. [lj. 


The ratio between non-perturbative constants C^ c /Cx^ N ^ for the model. The meaning of the symbols is the same 
Some data have been horizontally shifted to render the figure clearer. 
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FIG. 3. The PT ratio for the model. Empty circles (squares, diamonds, triangles) stand for the 1 -loop (2-loop, 3-loop, 
4-loop) results in the standard scheme; filled circles (squares, diamonds, triangles) stand for the 1 -loop (2 loop, 3-loop, 4 loop) 
results in the energy scheme. Some data have been horizontally shifted to render the figure clearer. 



FIG. 4. The KT ratio for the model. 
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TABLE I. Results of Monte Carlo data for the 0(3) nonlinear cr-model with standard action. 


p 

L 

stat 

X 


E 

1.50 

100 

5 x 10° 

176.20(0.16) 

10.99(1) 

0.6015813(19) 

1.60 

150 

5 x 10 6 

446.91(0.42) 

18.89(2) 

0.6357033(11) 

1.70 

260 

5 x 10® 

1264.4(1.00) 

34.36(3) 

0.6642223(6) 

1.75 

340 

5 x 10 6 

2189.5(2.2) 

46.97(5) 

0.6766299(4) 

1.80 

450 

5 x 10 6 

3827.4(3.9) 

64.43(7) 

0.6879333(3) 

1.85 

640 

5 x 10 6 

6717.6(7.1) 

88.53(10) 

0.6983241(2) 

1.90 

860 

5 x 10 6 

11850.(13.) 

121.70(10) 

0.7079167(1) 
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